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Abstract 

The breaking of conventional linear ^-factorization for hard processes in a nuclear envi- 
ronment is by now well established. Here we report a detailed derivation of the nonlinear 
/c^-factorization relations for the production of quark-gluon dijets. This process is of direct 
relevance to dijets in the proton hemisphere of proton-nucleus collisions at energies of the Rel- 
ativistic Heavy Ion Collider (RHIC). The major technical problem is a consistent description 
of the non-Abelian intranuclear evolution of multiparton systems of color dipoles. Following 
the technique developed in our early work [ N.N. Nikolaev, W. Schafer, B.G. Zakharov and 
V.R. Zoller, J. Exp. Theor. Phys. 97 (2003) 441], we reduce the description of the intranuclear 
evolution of the qggq state to the system of three coupled-channel equations in the space of 
color singlet 4-parton states |33), 1 66) and 1 1515) (and their large- N c generalizations). At 
large number of colors N c , the eigenstate ( 1 66) — |1515))/\/2 decouples from the initial state 
1 33). The resulting nuclear distortions of the dijet spectrum exhibit much similarity to those 
found earlier for forward dijets in Deep Inelastic Scattering (DIS). Still there are certain distinc- 
tions regarding the contribution from color-triplet qg final states and from coherent diffraction 
excitation of dijets. To the large-iV c approximation, we identify four universality classes of 
nonlinear fcj_-factorization for hard dijet production. 



*Electronic address: N.N. Nikolaev@fz-juelich.de 
^Electronic address: Wo.Schaefer@fz-juelich.de 
* Electronic address: B.Zakharov@fz-juelich.de 
^Electronic address: zoller@her on.itep.ru| 



1 



I. INTRODUCTION 



According to the conventional perturbative QCD (pQCD) factorization theorems the 
hard scattering cross sections are linear functionals (convolutions) of the appropriate 
parton densities in the projectile and target l|. An implicit assumption behind these 
theorems is that the parton densities in the beam and target are low and the relevant 
partial wave amplitudes are small, so that the unitarity constraints can be ignored. In 
the case of hard processes in a nuclear environment, the properly defined partial wave 
amplitudes become proportional to the nuclear thickness and, for a sufficiently heavy 
nucleus, overshoot the s-channel unitarity bound. The unitarization makes the nuclear 
partial waves a highly nonlinear functional of the free nucleon amplitudes. Alternatively, 
in the pQCD language, the unitarity constraints bring in a new dimensional scale into the 
problem - the so-called saturation scale. Important implication of the nonlinear unitarity 
relation between the free-nucleon and nuclear partial waves is that the properly defined 
density of gluons in a nucleus becomes a nonlinear functional of the gluon density in 
a free nucleon; the first discussions of the fusion of partons in deep inelastic scattering 

Q 

(DIS) off a nucleus go back to 1975 [2j|. 

The emergence of a new large scale and the ensuing nonlinearity call for a revision 
of the pQCD factorization for hard processes in a nuclear environment. A consistent 
analysis of forward hard dijet production in DIS off nuclei revealed a striking breaking of 
linear fcj_-factorization P, 4| confirmed later on in the related analysis of single-jet spectra 
in hadron-nucleus collisions 0,0]. Namely, following the pQCD treatment of diffractive 
dijet production 0,0], one can define the collective nuclear unintegrated gluon density 
such that it satisfies the s-channel unitarity constraints and such that the familiar linear 
fc_i_-factorization (see e.g. the recent reviews [2J) would hold for the nuclear structure 
function F 2 a(x, Q 2 ) and forward single-quark spectrum in DIS off nuclei because of their 
special Abelian features. However, the dijet spectra in DIS and single-jet spectra in 
hadron-nucleus collisions prove to be a highly nonlinear functionals of the collective 
nuclear gluon density. Furthermore, the pattern of nonlinearity for single-jet spectra 
was shown to depend strongly on the relevant partonic subprocess 0. Our conclusions 
on the breaking of linear fc^-factorization for hard scattering off nuclei were recently 
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taken over by other authors 
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In this communication we extend the analysis |4j,|£|,|13j of the excitation of heavy flavor 
and leading quark dijets in DIS, Y9n QQ, where stands for the gluon exchanged 
with the nucleon, to the excitation of quark-gluon dijets (pQCD Bremsstrahlung tagged 
by a scattered quark) in the pQCD subprocess q*gN Q9 off free nucleons and its gen- 
eralization to heavy nuclear targets. In the latter case multiple gluon exchanges between 
the involved partons and a nucleus are enhanced by a large nuclear radius. The issues 
are (i) to which extent such multiple gluon exchanges can be described in terms of the 
unintegrated collective nuclear gluon density and (ii) whether the nuclear factorization 
for quark-gluon dijets in qA collisions is similar to that for the quark- ant iquark dijets 
in DIS, i.e, in collisions. To a certain extent, our answer is in the affirmative - the 
nonlinear /c^-factorization properties for two processes exhibit much similarity. Still, the 
two cases differ substantially. For instance, the production of coherent diffractive dijets 
makes about 50% of the total cross section in DIS but becomes marginal in qA collisions. 
Furthermore, the contributions from quark-gluon dijets in different color multiplets have 
a very distinct nonlinear /^-factorization properties. Also the effects of the initial state 
interaction change substantially from the color-singlet 7* in DIS to the color-triplet quark 
in qA collisions. On the other hand, the unifying aspect is a treatment of the excitation 
of final-state color dipoles in the higher color multiplets - color-octet in DIS and sextet 
and 15-plet in qA collisions. 



for the inclusive dijet 
14 , [3| and allows to 



The starting point of our analysis is the master formula 
spectrum. It is derived based on the technique developed in 
calculate the dijet spectrum in terms of the S-matrices for interaction with the target 
nucleon or nucleus of the color-singlet n-parton states, n = 2, 3, 4. Within this technique, 
one deals with infrared-safe quantities despite the fact that the incident parton - the 
quark q* - is carrying a net color charge. The calculation of the two-parton and three- 
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parton S-matrices is the single-channel problem with the known solution 
stumbling block is the calculation of the 4-body S-matrix. In the case of the quark-gluon 
dijets it describes the non-Abelian intranuclear evolution of the color-singlet qgqg system 
of dipoles. It can be reduced to a 3 x 3 coupled-channel problem. In our earlier work j^J 
we published a full solution of the related two-channel problem for the qqqq system which 



emerges in the description of quark-antiquark dijets in DIS. Here we report the solution 
for the qgqg system of dipoles which has some new features compared to the qqqq state 
in DIS 1 . We go to fine details of this derivation - specifically, the diagonalization of the 
coupled-channel S-matrix and to the formulation of explicit nonlinear fc^-factorization 
formulas for the dijet spectrum - for several reasons. First, the production of quark- 
gluon dijets without the soft gluon approximation has not been treated before. Second, 
regarding the color properties of the incident and final states, it is a process of sufficient 
generality to set a basis for the description of other pQCD processes. In conjunction 
with our earlier results, it allows to formulate four universality classes of nonlinear k±- 
factorization. Third, recently the formal representation for the dijet cross section similar 
to our master formula has been discussed by several groups [r], [ill Q]; but these 
works stopped short of the diagonalization of their counterpart of our 4-body S-matrix. 
Correspondingly, they do not contain explicit nonlinear fc^-factorization formulas. 

A very rich pattern of the pro cess- dependent nonlinear /cj_-factorization emerges from 
the studies presented here and reported in J, 0, 13, [r|. For instance, it becomes 
increasingly clear that a heavy nucleus cannot be described in terms of a universal 
collective glue, rather the nuclear glue must be described by the density matrix in the 
color space. Furthermore, the collective glue defined for the slice of a nucleus rather than 
the whole nucleus is an integral part of the description of excitation of color dipoles in 
higher color representations. The linear /c^-factorization for the single-quark jets in DIS 
found in our earlier study [4] is an exception due to the Abelian incident parton - the 
photon. 

From the point of view of practical applications, the discussed quark-gluon dijets are 
of direct relevance to the large (pseudo) rapidity region of proton-proton and proton- 
nucleus collisions at RHIC (for the discussion of the possible upgrade of detectors at 
RHIC II for the improved coverage of the proton fragmentation region see Our 
treatment is applicable when the beam and final state partons interact coherently over 
the whole longitudinal extension of the nucleus, which at RHIC amounts to the proton 



1 A brief discussion of the main results has been reported elsewhere [lfij . 
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fragmentation region of 



(Q*) 2 + M? 1 „ l7 o 

„ ^ X < x A = — « 0.1A- 1/3 , 



2mE g * ~ ~ 2R,ATn, 
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where i?^ is the radius of the target nucleus of mass number A, (Q*) and E q * are 
the virtuality and energy of the beam quark q* in the target rest frame, M± is the 
transverse mass of the dijet and m p is the proton mass (U f° r the related color 
dipole phenomenology of the experimental data on nuclear shadowing see |12|). 

The presentation of the main material is organized as follows. The master formula 
for the dijet spectrum is presented in Sec. II. The two-body density matrix - the Fourier 
transform of which gives the dijet spectrum - contains the S-matrices for the interaction 
of two-, three- and four-parton color-singlet systems of dipoles with the target. Based on 
the technique developed in 1 151, in Sec. Ill we report single-channel S-matrices in terms of 

nn 

the quark- ant iquark and quark-antiquark-gluon color-dipole cross sections [lo,l2I|- Sec. 
IV contains all the technicalities of the derivation of the coupled-channel S-matrix for the 
qgqg state: the decomposition into color multiplets; projection onto the final states; the 
color-flow diagram technique for the calculation of the 3x3 matrix of color-dipole cross 
sections; the derivation of the relevant Casimir operators; the explicit diagonalization of 
the S-matrix at large number of colors N c and the Sylvester expansion. The quark-gluon 
dijet spectrum for the free-nucleon target is derived in Sec. V. Here we also comment on a 
direct relationship between the dijet and single-jet spectra for the free-nucleon reactions 
described by the single-gluon exchange in the t-channel. The principal result of this study 
- the nonlinear fc_|_-factorization for the dijet spectrum produced off nuclear targets - is 
reported in Section VI. Here we compare the pattern of nonlinear /c^-factorization for 
quark-gluon dijets in qA collisions to that for the quark-antiquark dijets in DIS and gA 
collisions and identify four universality classes depending on the color representation of 
the incident parton and final-state dijet. In Section VII we apply our results to the 
nuclear broadening of the dijet acoplanarity distribution. In Sec. VIII we comment 
on a limiting case when the quark-gluon dijets merge to one jet. Such monojets can be 
identified with the fragmentation of the quark jet formed by the quasielastically scattered 
incident quark. The separation into the dijet and monojet final states changes with the 
mass number of the target nucleus and the centrality of the collisions. We also comment 
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FIG. 1: The rapidity structure of the radiation of gluons by quarks q — > qg in pA collisions. 

on the possible nuclear modification of the fragmentation function. In the Conclusions 
section we summarize our main results. 

II. THE MASTER FORMULA FOR QUARK-GLUON DUET PRODUCTION 
OFF FREE NUCLEONS AND NUCLEI 

A. Kinematics and nuclear coherency 

To the lowest order in pQCD the underlying subprocess for quark-gluon dijet produc- 
tion in the proton fragmentation region of proton-nucleus collisions is a collision of the 
quark q* from the proton with the gluon from the target, 

q*9N -> qg ■ 

It is a pQCD Bremsstrahlung tagged by a scattered quark. We don't restrict ourselves to 
soft gluons. In the case of a nuclear target one has to deal with multiple gluon exchanges 
which are enhanced by a large thickness of the target nucleus. 

From the laboratory, i.e., the nucleus rest frame, standpoint it can be viewed as an 
excitation of the perturbative \qg) Fock state of the physical projectile \q*) by one-gluon 
exchange with the target nucleon or multiple gluon exchanges with the target nucleus. 
Here the collective nuclear effects develop if the coherency over the thickness of the 
nucleus holds for the qg Fock states, i.e., if the coherence length is larger than the 
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diameter of the nucleus, 



lc = ~ypt — > 2 -^< 

Mi xttin 



where 

M 2 ± = P 1 + P R (3) 

Zq Z g 

is the transverse mass squared of the qg state, p qg and z q ^ g are the transverse momenta 
and fractions of the the incident quarks momentum carried by the quark and gluon, 
respectively (z q + z g = 1). In the antilaboratory (Breit) frame, the partons with the mo- 
mentum xpn have the longitudinal localization of the order of their Compton wavelength 
A = 1/xpN, where Pat is the momentum per nucleon. The coherency over the thickness 
of the nucleus in the target rest frame is equivalent to the spatial overlap of parton fields 
of different nucleons at the same impact parameter in the Lorentz-contracted ultrarel- 
ativistic nucleus. In the overlap regime one would think of the fusion of partons form 
different nucleons and collective nuclear parton densities [2]. The overlap takes place if 
A exceeds the Lorentz-contracted thickness of the ultrarelativistic nucleus, 

\=^->2R A -?», (4) 
xp N p N 

which is identical to the condition (J2J). 

Qualitatively, the both descriptions of collective nuclear effects are equivalent to each 



1 devel- 
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other. Quantitatively, the laboratory frame approach takes advantage of the we 
oped multiple-scattering theory of interactions of color dipoles with nuclei 4j,[l8, 
From the practical point of view, the coherency condition x < xa restricts collective 
effects in hard processes at RHIC to the proton fragmentation region. The target 
frame rapidity structure of the considered q* — > qg excitation is shown in Fig. ^ The 
(pseudo)rapidities of the final state partons must satisfy r] qtg > t]a = logl/x/i- The 
rapidity separation of the quark and gluon hard jets, 

1 — z 

Ar]qg = log a - , (5) 

Z 9 

is considered to be finite. Both jets are supposed to be separated by a large rapidity 
from other jets at mid-rapidity or in the target nucleus hemisphere; the gaps between 
all jets, beam spectators and target debris are filled by soft hadrons from an underlying 
event . 
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B. Master formula for excitation of quark-gluon dijets 



In the nucleus rest frame, relativistic partons q*,q and g, propagate along straight- 
line, fixed-impact-parameter, trajectories. To the lowest order in pQCD the Fock state 
expansion for the physical state \q*) p h ys reads 

\q*)phys = \q*)o + y(zg,r)\qg} , (6) 

where ^{z g , r) is the probability amplitude to find the qg system with the separation r 
in the two-dimensional impact parameter space, the subscript "0" refers to bare partons. 
The perturbative coupling of the q* — > qg transition is reabsorbed into the lightcone wave 
function ty(z g ,r). We also omitted a wave function renormalization factor, which is of 
no relevance for the inelastic excitation to the perturbative order discussed here. The 
explicit expression for ^(z g , r) in terms of the quark-splitting function will be presented 
below. The wave function depends on the virtuality of the incident q*, which equals 
(Q*) 2 = (p*) 2 , where p* is the transverse momentum of q* in the incident proton (Fig. [I}. 
For the sake of simplicity we take the collision axis along the momentum of the incident 
quark q*, the transformation between the transverse momenta in the g*-target and p- 
target reference frames is trivial 

If b is the impact parameter of the projectile q*, then 

b q = b- z g r } b g = b + z q r. (7) 

By the conservation of impact parameters, the action of the S- matrix on \a) p h ys takes a 
simple form 

S\q*) P hys = S q {b)\q*) + S q {b q )Sg(b g )ty{z,r)\qg)o 

= S q (b)\q*) phys + [S q (b q )S g (b g ) - S q (b)}*(z g ,r)\qg) . (8) 

In the last line we explicitly decomposed the final state into the (quasi)elastically scat- 
tered \q*)ph ys and the excited state \qg)o- The two terms in the latter describe a scatter- 
ing on the target of the qg system formed way in front of the target and the transition 
q* — > qg after the interaction of the state |g*)o with the target, as illustrated in Fig. 
121 The contribution from transitions q* — > qg inside the target nucleus vanishes in the 
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a) b) c) 

FIG. 2: Typical contribution to the excitation amplitude for q*A — ► qgX, with multiple color 
excitations of the nucleus. The amplitude receives contributions from processes that involve 
interactions with the nucleus after and before the virtual decay which interfere destructively. 

high-energy limit of a; < xa 2 ■ We recall, that the s-channel helicity of all partons is 
conserved. 

The probability amplitude for the two-jet spectrum is given by the Fourier transform 

J d 2 b q d 2 b g exp{-t(p q b q+Pg b g )][S q (b q )S g (b g ) - S 9 (b))*(z g ,r) (9) 

The differential cross section is proportional to the modulus squared of (j§J, 

/ d%d%eMi(p q K + P 9 b' g )][Sl(b' q )Sl(b' g ) - S q (b'W(z g ,r') 
x J d 2 b q d 2 b g exp[-t(p q b q+Pg b g )}[S q (b q )S g (b g ) - S q (b)Mz b ,r) . (10) 

The crucial point is that the hermitian conjugate can be viewed as the S-matrix for 
an antiparton jj, 14, Consequently, the four terms in the product 

[s,(&ys,(6y - s q (b')ns q (b q )s g (b g ) - s q (b)\ 

admit a simple interpretation: 

S«> (l/,6) = SUb')S q (b) (11) 



2 In terms of the lightcone approach to the QCD Landau-Pomeranchuk-Migdal effect, this corresponds 
to the thin-target limit |22j. 
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can be viewed as an S-matrix for elastic scattering on a target of the q*q* state in which 
the antiparton q* propagates at the impact parameter b' . The averaging over the color 
states of the beam parton q* amounts to the dipole q*q* being in the color singlet state. 
Similarly, 

Sf q9 (b',b q ,b g ) = Sl(b')S g (b q )S g (b g ), 

s %A b ^K) = si(b' q p g (b' g )s q (b) 

Sf 9 > 9q K K> K b g ) = Sl(b' q )Sl(b' g )S g (b g )S q (b q ) . (12) 

describe elastic scattering on a target of the overall color-singlet qqg and qggq states, 
respectively. This is shown schematically in Fig. EJ Here we suppressed the matrix 
elements of over the target nucleon, for details of the derivation based on the closure 
relation, see 4J]. Specifically, in the calculation of the inclusive cross sections one averages 
over the color states of the beam parton q*, sums over color states X of final state partons 
q,g, takes the matrix products of and S with respect to the relevant color indices 
entering and sums over all nuclear final states applying the closure relation. The 
technicalities of the derivation of will be presented below, here we cite the master 
formula for the dijet cross section, which is the Fourier transform of the two-body density 
matrix: 

dzd 2 p q d?p g (2n) 4 J 1 9 b c 
x exp[-ip q (b q - b' q ) - t Pg (b g - b' g ))*(z g , b q - b g )V*(z g , b' q - b' g ) (13) 

^(X\^ g ) qg (b' q , b' g , b q , b g ) + S%{b', b) - Sf g ] q ,(b, b' q , b' g ) - Sg g (b', b q , 6 ff )||m) 

Hereafter, we describe the final state dijet in terms of the gluon jet momentum, 
p = p g , z = z g , and the decorrelation (acoplanarity) momentum A = p q + p g - We also 
introduce 

s = b q -b' q , (14) 
in terms of which b g — b' g = s + r — r' and 

exp[-ip q (b q - b' q ) - ip g (b g - b' g )] = exp[-iAs - ipr + ipr') , (15) 

so that the dipole parameter s is conjugate to the acoplanarity momentum A. 

10 



4 



+ 






q (4) + q (2) 
Oqgqg Oq* q " 






(3) 



S V _ 
a 



qgq 



q(3) 

^qgq' 



FIG. 3: The S-matrix structure of the two-body density matrix for excitation q* 
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III. CALCULATION OF THE 2-PARTON AND 3-PARTON S-MATRICES 

A. The quark-nucleon S-matrix and the /c^-factorization representations for 
the color dipole cross section 



In order to set the formalism, we start with the S-matrix representation for the cross 
section of interaction of the triplet-antitriplet color dipole qq with the free-nucleon target 
Pi . To the two-gluon exchange approximation, the S-matrices of the quark-nucleon and 
antiquark-nucleon interaction equal, respectively, 

S(6 fl ) = 1 + iT a V aX (b q ) - l -T a T a x \b q ) , 

&(b 9 ) = 1 - iT a V aX (b,) - l -T a T a X 2 (b q ) , (16) 

were T a V a x{b) is the eikonal for the quark-nucleon gluon exchange. The vertex V a 
for excitation of the nucleon g a N — * N* into color octet state is so normalized that 
after application of closure over the final state excitations N* the vertex g a g b NN equals 
(N\V^Vb\N) = S a b. The second order terms in ([To]) do already use this normalization. 
The S-matrix of the (gg)-nucleon interaction equals 

c(%, , n (N\TT[S(b q )S\b,)]\N) 



[b q , bg 



(17) 



(N\1\N)Tt1 

A graphical rule for the calculation of the color traces entering (fTTj) is shown in Fig. 0] 
such color flow diagrams will extensively be used in the subsequent calculation of S^ 4 -*. 
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S m (a-a)= 



S q (a) 



S*(a) 




I 



FIG. 4: The color-flow diagram for the S-matrix for the interaction of the color-single qq dipole 
with the nucleon; a and a are the impact parameters of the quark and antiquark, respectively. 



(2) 

The corresponding profile function is T 2 (b q , bg) = 1 — S qq (b q ,b q ). The dipole cross 
section for interaction of the color-singlet qq dipole r = b q — b q with the free nucleon is 
obtained upon the integration over the overall impact parameter 

a{r) = 2 J d 2 b q T 2 (b q , b q -r)=C F J d%[ X (b q ) - X (b q - r)] 2 , (18) 

where Cp = (N 2 — 1)/2N C is the quark Casimir operator. It sums a contribution from 
the four Feynman diagrams of Fig. El and is related to the gluon density in the target 
by the /c^-factorization formula |2JJ, |2^ 

a(x, r) = J d 2 nf(x, n)[l — exp(inr)] , (19) 

where 

ft \ 47ra 5 (r) 1 2 

f(x,K) = —j^ — -T{x,k ) (20) 

and 



dG(x, k 2 
d log K? 

is the unintegrated gluon density in the target nucleon. Hereafter, unless it may cause a 
confusion, we suppress the variable x in the gluon densities and dipole cross sections. The 
leading Log^ evolution of the dipole cross section is governed by the color-dipole BFKL 
evolution [20, l^^j , the same evolution for the unintegrated gluon density is governed by 
the familiar momentum-space BFKL equation 
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FIG. 5: The four Feynman diagrams for the quark-antiquark dipole-nucleon interaction by the 
two-gluon pomeron exchange in the t-channel. 



The S-matrix for coherent interaction of the color dipole with the nuclear target is 
given by the Glauber-Gribov formula j3, E 



where 



S[b,a(r)} = exp[~a(r)T{b)] 



/oo 
dr z n A {b, r z 
-oo 



(22) 



(23) 



is the optical thickness of the nucleus. The nuclear density riA{b,r z ) is normalized 
according to J d 3 fn^(b,r z ) = J d 2 bT(b) = A, where A is the nuclear mass number. 



In the specific case of S-J t (b', b) the color dipole equals 



99 



Tqq = b — b' = s + zr — zr 



(24) 



and S^j? ,(i/, b) entering Eq. (fT4*|) will be given by the Glauber-Gribov formula 



S%*(b',b)=S[b,*(s + zr-zr')) 



(25) 



B. The S-matrix for the color-singlet qqg state 



Here quark and gluon couple to the co 



or triplet. The dipole cross section for the 

n 

color singlet qqg state has been derived in |20f] , the S-matrix derivation with the quark- 



antiquark basis description of the gluon is found in Appendix A of ref. [6]. For the 
generic 3-body state shown in Fig. |U1 it equals 

C 

°3(r qq -, r gq ) = 2£-[cr(r gq ) + <r(r gq ) - a(r qq )] + a{r qq ) , (26) 
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a) b) 

FIG. 6: The color dipole structure of (a) the generic quark-antiquark-gluon system of dipoles 
and (b) of the q*qg system which emerges in the S-matrix structure of the two-body density 
matrix for excitation q* — ► qg. 



where r gq = r gq + r qq . The configuration of color dipoles for the case of our interest is 
shown in Fig. El (see the related derivation in jl5|). For the q*qg state the relevant dipole 
sizes in equal 



r qq 


= b q 


— b' = s — zr , 




r gi 


= b 9 


-b q = r, 






= b g 


— b' = s + r — zr' , 


(27) 



whereas for the q*qg' state we have 
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b - b' q = s + zr, , 

b' 9 -b' q = r\ 

b g — b = s + zr — 



(28) 



so that 



Q 'qg 



&q*qg' 



C A 
2C F 
Ca 
2C f 



[cr(r) + cr(s + r — zr') — cr(s — zr')] + cr(s — zr) , 
-[a(-r') + cr(s — r' + zr) — a(s + zr)] + cr(s + zr) 



(29) 



The overall color-singlet qqg state has a unique color structure and its elastic scattering 
on a nucleus is a single-channel problem. Consequently, the nuclear S-matrix is given by 



14 



the single-channel Glauber-Gribov formula |26l . |2 

Sf g > q *(b,b' g ,b g ) = S^^W 



S§L(&'AA) = S[b,a rw ]. (30) 



IV. COUPLED-CHANNEL S-MATRIX FOR THE 4-PARTON STATE 
A. The basis of color-singlet (qqgg 1 ) states 

The 4-parton S-matrix describes transitions between color-singlet (qqgg') states. It 
is convenient to decompose the the \qg) state into the |3), |6) and |15) states and 
their SU(N C ) generalizations (our reference to the triplet, sextet and 15-plet states at 
arbitrary N c should not cause any confusion). Then the basis of color-singlet states 
\qqgg') will consist of the |33), 1 66) and 1 15 15) systems of color dipoles and the in- 
tranuclear evolution in the elastic scattering of the 4-parton state off the nucleus is the 
three-channel problem. The evolution starts from the 1 33) state what is evident from 
Fig. El Technically, once the 3x3 matrix £ of 4-body dipole cross sections is known, 
the corresponding nuclear S-matrix will be given by the Glauber-Gribov formula 

S ( H, gq (b' q ,b' g ,b q ,b g ) = S[b,t]. (31) 

Our immediate task is a calculation of the coupled-channel operator £. 
We chose a description of the gluon in the quark- ant iquark basis: 

g[ = a i a k - ^[aa)8i. (32) 

In the calculation of the S-matrices both the quark a and the antiquark a must be 
considered as propagating at the same impact parameter. The generic quark-gluon state 
is described by a tensor 

4/ = 9k c i = a l a k ci - -^-(aa)a5l . (33) 

There is a unique color-triplet quark-gluon state (the normalization of the states will be 
defined at the level of the |33), 1 66) and |15 15) systems of color dipoles) 

t k = [dc)a k - -^-{aa)c k (34) 
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The sextet state is described by the traceless tensor antisymmetric in (k,l), 

A li = a l {a k Ci - aic k ) + — — - [(ac)a f - {aa)a}5' k - — — - [{ac)a k - {aa)c k }5\ , (35) 

while the 15-plet is described by the traceless symmetric tensor 

S kl = a l (a k ci + aic k ) - ^— - [(ac)ai + (aa)ci]5 l k - — — ^ [{ac)a k + (aa)c k ]5\ . (36) 

The quark, antiquark and two gluons in the color-singlet (qqgg f ) system of dipoles 
all propagate at different impact parameters. To avoid a confusion, the gluon in the 
complex conjugated state will be described by the tensor 



and the antitriplet state is 



(g% = b%-±(bb)5i 



t k = (cb)b k - ^{bb)c k . 



(37) 



(38) 



The overall color-singlet |33), 1 66) and |15 15) states will be decomposed into six 6-body 
color-singlet states. The corresponding 6-body vertices (projection operators) equal 



V\ = (ab)(ba)(cc), Vi = (ab)(bc)(ca), V3 = (aa)(bc)(cb), 
V4 = (ac)(bb)(ca), V5 = (ac)(ba)(cb), Vq = (aa)(bb)(cc). 



(39) 



some of which are pictorially represented in Fig. [7| For instance, the normalized color- 
singlet triplet-antitriplet state will be 



|33) 

Similarly, one finds 



(JV? - 1) 



(40) 



|66)= V 1 -V 2 + 



N c -1 



(V 3 + v 4 -v 5 - v 6 \ 



2N c (N c + l)(N c -2) 



(41) 



1 15 15) 



V 1 + V 2 



1 



iV c + l 



(V 3 + v 4 + v 5 + v 6 



1 



2N c (N c -l)(N c + 2) 



(42) 



These states are normalized to unity, (33 1 33) = (66|66) = (15 15|15 15) = 1, the normal- 
ization coefficients are readily derived using the color-flow diagram technique described 
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a 
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a 



V^-{ab){ba){cc) V 2 = (ab)(ba)(cc) V 5 = (ac)(ba)(cb) 

FIG. 7: Examples of the 6-body vertices (projection operators) which emerge in expansions of 
the qgqg states in the quark-antiquark basis. 



in Sec. IV-C below. The diagonal and off-diagonal matrix elements of the 4-body cross 
section operator in the basis of 1 33) , 1 66) and |15 15) of color dipole states will be decom- 
posed in terms of the matrix elements 

<r* = (VMVk) (43) 

with the coefficients which are readily read from the expansions (I4U|) - (|42I) . 

Note, that each of the is a matrix element between the overall color-singlet 6-body 
configurations composed of the three color-singlet quark-antiquark pairs. As such all of 
them are infrared-safe quantities. 



B. Projection onto the final states 

In the case of the inclusive dijet spectrum with summation over all colors of final state 
quarks and gluons the projection onto the final state is of the form (see the discussion 
in Q) 



J2( X \ = E^im(R)(RR\ 

X R 



■V,(33| + J\n c (N c + 1)(N C - 2)(66| + J\n c (N c - 1)(N C + 2)(15 15| , (44) 
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where dim(i?) is the dimension of the corresponding representation. The averaging over 
the colors of the initial quark q* amounts to taking 

1 



\in) = |33) 



(45) 



'dim^) 

Then the calculation of the inclusive dijet cross section requires the evaluation of the 
combination of matrix elements 



J2(X\S[b,E]\in) = (33|S[6,S]|33) + 
x 



+ 



dim(6) 



(66|S[6, E]|33) 



dim(15) (15T5|S[6,S]|33) 



(46) 



\ dim(3) x 1 L ' Jl ' \\ dim (3) 

Besides the inclusive spectrum one can readily consider the excitation of quark-gluon 
dijets in specific color representations. We reiterate that they also will be infrared-safe 
observables. 



C. Color-flow diagrams 

The calculation of the matrix elements ()43j) is greatly simplified by the technique 
of color-flow diagrams. Each matrix element (J43j) corresponds to a certain color flow 
diagram. Altogether there are 21 different color flow diagrams, the three selected cases 
are shown in Fig. |HJ The number of closed loops varies from three to one. In the 
calculation of the S-matrix elements, 

S ik = (Vi\S\V k ) , (47) 

each horizontal quark line is multiplied by the quark S-matrix S(bj) taken at the appro- 
priate impact parameter bj, while the antiquark line is multiplied by S^(bj). The trace 
of the product of S-matrices is calculated for each closed loop. 

The first application of the color-flow diagrams is to the derivation of the normaliza- 
tion factors in expansions (}4~Tj) . They are obtained by assigning the factor N c to each 
and every loop. 

Now we present the results for the three matrix elements shown in Fig. |H| For the 
sake of brevity the impact parameters of quarks and antiquarks will be denoted by their 
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^11 w 12 w 25 

FIG. 8: Examples of color flow diagrams for the calculation of the components of the 6-body 
dipole cross sections. The horizontal quark lines are multiplied by the quark S-matrix S(foj) 
taken at the appropriate impact parameter, while each horizontal antiquark line is multiplied 
by S^(bj), the trace is taken for each closed loop. 



symbols. One readily finds 
Su = (K|S|K) 



Tr 



S(a)S t (6)l Tr[S(6)S t 



Tr 



S(c)St(c) 



at3[i _ r( fl - 5)][i - r(6 - a)][i - r( c - c)] 



(4£ 



The multibody S-matrices must be evaluated up to the terms quadratic in the QCD 
eikonal, i.e., to the terms linear in the triplet-antitritplet color-dipole profile function T, 
and the corresponding matrix element of the dipole cross section equals 



a n = (Vtla^Vt) = N^[a(a - b) + a{b - a) + a(c - c)] 
= N*[2a{a - b) + cr(c - c)) . 



(49) 



Each quark-antiquark loop gives the corresponding dipole cross section, times N c to the 
power equal to the number of loops. Here we took into account that the quark a and 
antiquark a, and b and b as well, propagate pairwise at equal impact parameters. 
The case of <j\2 is a bit more complicated. Here Si 2 is a product of two traces: 



S12 = (^ 1 |S|^ 1 ) = Tr[S(fe)S t 



Tr 



S^S^StcOS^c) 



N e [l - r(6 - a)]Tr[S(a)S t (6)S(c)S t (c 



(50) 
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b) 



FIG. 9: Examples of interaction with the target nucleon of the (a) quark-antiquark and diquark 
dipole in the bcca state 



The latter trace Tr S(a)S^(6)S(c)S^(c) was already encountered in our derivation of 
the 4-parton S-matrix for the production of dijets in DIS [4]. The corresponding color- 
flow diagram is shown in Fig. El Here one needs to sum the contributions to the bcca 
scattering amplitude from the exchange by the 2-gluon pomerons in the t-channel. The 
familiar diagram of Fig. Ek gives a contribution — %(c)x(c)Tr(T a T a ) . The new case is 
when the two gluons are attached to the diquark ac as shown in Fig. Eb- Straightforward 
color algebra shows that the corresponding contribution to the profile function equals 
x(a)x( c )Tr(T a T a ) . This gives rise to a simple rule: each quark-antiquark pair, ab, ac, cb 
and cc, contributes the corresponding triplet-antitriplet dipole cross section, whereas the 
diquark ac and the anti-diquark be contribute the triplet-antitriplet dipole cross section 
taken with the negative sign. The color traces give a factor N c per each loop, one of 
these factors has already been put in evidence in Eq. (JBTlj) . The final result is 

tTi2 = (ViH^i) = N 2 c [o(b-a) + o(a-b) 

— a(a — c) + a (a — c) + er(c — b) — a(b — c)) + cr(c — c)] . (51) 

Application of the same technique to cr 2 5 gives 



'25 



(K|S|Vi) = Tr[S(a)S t (a)S t (c)S(&)S t (a)S(c)S t (&)] 



(52) 
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with the cross section 

025 = (V 2 \(T A \V 5 ) = N c [a(a - c) - a(a - b) 

+a(a — a) — a(a — c) + a(a — b) + a(b — c) — a(a — c) 

+a(c - c) - a(c -b) + a(b -a)- a(b - c) + a(b - b) 

+a(c -a) - a(a -b) + a(c - b)] = N c a(c - c) (53) 

Here we used the obvious properties a (a — a) — a(b — b) = and cancellations due to 
equalities of the form a(c — a) = a(c — a). For the sake of completeness, we cite all the 
remaining a^: 



ova 


= AT c a(c- 


c), 




(Ji4 


= iV c cr(c- 


c), 




Cl5 


= N 2 c [2a(a 


- b) + <j(c - 


-c) + 




+ cr(a — c) 


+ <r(b - c) - 


- (7(6 - c 


0"16 


= iV c V(c - 


■c), 




0"22 


= N?[cr(a- 


-b) + a(b - 


c) + <r(a 


0"23 


= N 2 c [a(b- 


- c) + a(b - 


c)], 


024 


= N 2 c [a{a- 


- c) + a(a - 


c)}, 


0"26 


= N c a(c- 


c), 




0"33 


= N 3 c [a(b- 


- c) + a(b - 


c)}, 


0"34 


= N c a(c- 


c), 




0"35 


= N?[a(b- 


- c) + a(b - 


c)], 


^36 


= iV c V(c - 


■c), 




<7 44 


= Nl[a{a- 


- c) + a(a - 


c)}, 


045 


= N 2 c [a(a- 


- c) + a(a - 


■ 5)] , 


046 


= N 2 a(c ~ 


■c), 




055 


= N 3 c [a(a- 


-b) + a(a - 


- c) + tr(6 


056 


= N c a(c- 


c), 




066 


= N?a(c - 


c). 
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D. The 3x3 matrix of 4-parton dipole cross sections S. 



A simple algebra gives the following 3x3 matrix £ of 4-body dipole cross sections 
(we go back to the dipole parameters defined in Sect. 2): 

C A 



S n = (33|a 4 |33) 
1 



N 2 -l 



a s - 



2C F 
C A 



<j(s — r' + r) + <r(r) + er(r')] 

1 n. 



2C F N 2 -l ' 



where 



f2 = cr(s — r') + cr(s + r) — a(s — r' + r) — cr(s) . 



Similar calculation gives 



S 22 = (66|a 4 |66) 



3iV c + 1 1 



+ 



N c + 1 2 



2(i\T? - 1) 
+ J^TjHr) + a(r')] 



[cr(s — r' + r) + cr(s) 
a(s — r' + r) — <x(s)] 



2(iV c + 1) 



1 + 



1 



(n c - iy 



n 



S 33 = (1515|<7 4 |15l5) 



= ^ ' \ ■ [ a ( s ~r f + r) + <r(s)] 

N 2 + 1 

+ 2(A^ - 1) ^ [(7(S ~ r ' + r) ~ a(s)] 
- J^Hr) + a(r')] 



1 + 



1 



(n c + iy 



n. 



(55) 



(56) 



(57) 



(58) 



2(N C - 1) 

All the off-diagonal matrix elements for transition between different color representations 
are proportional to Q: 

N 2 



J 21 



(66 1 a 4 1 33) 



TV — 2 

n, 



(7V c -l)(iV c 2 -l)V2(iV c + l) 



(59) 



S 31 = (15 15|cr 4 |33> 



N 2 



N c + 2 



(N C + 1)(N 2 -1)^2(N C -1) 



n 



(60) 



22 




FIG. 10: The color- flow diagram representation of the quark Casimir operator Cf- 



S 32 = (1515|a 4 |66) 



1 

2 



-n. 



(61) 



Note, that the off-diagonal Q has precisely the same color-dipole structure as the off- 
diagonal o"i 8 which describes the excitation qq dipole from the color-singlet to color-octet 
state [4]. This off-diagonal matrix element vanishes if either r = of r' = 0, when the 
pointlike \qg) and \q'g') Fock states cannot be resolved. 



E. The pointlike triplet, sextet and 15-plet dipoles and Casimir operators 

In the limit of r = r' = 0, the 4-body states reduce to the pointlike triplet-antitriplet, 
sextet-antisextet and 15-15 dipoles. 



Indeed, in this limit 



as expected, while 



a(s) (62) 



3iV c + l 
S22 = cr(s), 

3N C — 1 

S33 = A^T a(s) • (63) 

The Feynman diagrams of Fig.^lmake it obvious that for partons in the representation R 
the dipole cross section must be proportional to the Casimir operator Cr. Consequently, 
the coefficients in (|63|) must equal the ratio Cr/Cf (a factor Cp for the triplet-antitriplet 
color dipole had been absorbed into the definition of cr(s), see Eq. Q18|) ). The derivation 
of Cr by the color-flow diagram technique proceeds as follows: 
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We recall that the calculation of Cp for the quark. 



_ Tr(r°r°) 

Trl ' V 7 



can be represented in terms of traces of color loop diagrams as shown in Fig. [TQj In 
order to avoid a confusion in the description of the conjugate states, it is convenient to 
represent the sextet qg state in terms of the three different quark fields, 

A li = a l {b k ci - b lCk ) + — l _ - [(ac)bi - (ab)ct]5 k - — l _ - [{ac)b k - {ab)c k \8\ . (65) 



One readily finds that 



AA oc (aa)(bb)(cc) — (aa)(bc)(cb) 

+ N l _ - (ac)(ba){cb) + — l _ - (aa)(bb)(cc) 

'ac)(bb){ca) — — (ab)(ba)(cc) . (66) 



N c -V /v /v ' N c -1 
In the quark representation the Casimir operator equals 

(T 6 + T c - T a ) 2 = 3C F + 2(T b T c ) - 2(T a T b ) - 2(T c T a ) . (67) 



The six color-flow diagrams generated by the expansion (J6TI|) are shown in Fig. El The 
straightforward calculation of the corresponding traces, putting the Tj on the relevant 
horizontal lines in the loops gives 

The similar expansion for the 15-plet state reads 

SS oc (aa)(bb)(cc) + (aa)(bc)(cb) 

(ac)(ba)(cb) + — -(aa)(bb)(cc) 



N c + V /v /v ' N c + 1 
+ ]^^(« c )(^)(^) + ^^(^)(H(cc) (69) 



and 



Cu-^C*. (70) 



This completes the check of the formulas (|63|) . 
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FIG. 11: The color-flow diagrams for the derivation of the Casimir operator Cf for sextet and 
15-plet qg states in the quark-antiquark representation. 



F. The N c — > — N c transformation between the sextet and 15-plet matrix ele- 
ments 

As a function of N c , the Casimir operators and matrix elements for transitions con- 
taining the sextet and 15-plet states satisfy a curious symmetry 

C 15 (N C ) = C 6 (-N c ), 
S 33 (iV c ) = a 22 (-N c ), 

S 13 (iV c ) = -a l2 (-N c ). (71) 

Evidently, the relative minus sign in the last line of (fTTj) is a matter of convention for the 
basis states. We do not offer any straightforward group-theoretic explanation for this 
transformation (see, however, a discussion of the correspondence between the symmetric 
and antisymmetric representations in Cvitanovic's lectures Q). 
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G. Large-iV c properties of £ 



The application of the above derived £ to the dijet production of the free-nucleon 

target is straightforward. In the case of the nuclear target one has to solve the secular 

equation for the eigenvalues and eigenstates of S. It is a cubic equation, can be solved 

in radicals and the corresponding eigenfunctions are directly calculable. The further 

fl 

application of the Sylvester expansion |4j to (|31|) is straightforward. Unfortunately, 
because of the radicals the relevant Fourier transforms in (pPfj) can only be performed 
numerically. Simple algebraic formulas for eigenvalues and analytic results for the dijet 
spectra are, however, obtained in the large-iV c approximation. The higher order terms 
of expansion in inverse powers of N c can also be presented in an analytic form ^j]. 
Note, that for large N c 

1 



£31 


— £21 


£32 


= -Q, 


2 ' 


S33 


= £ 2 2 



N c y/2 



2a(s 



n 



a s 



(72) 



which shows that one must first diagonalize the matrix £ in the |66), |15 15) sector. The 
two eigenvalues are 



J 2,3 



022 ± 



(73) 



and the corresponding eigenstates are 



1 



|2> = ^=(|66) 
|3> = ^(|66) 



]A 

|15 15>) 



3/2 



™, Vo 



|15 15>) 



3/2 • 



(74) 



In the basis of the states |1) = |33), |2) and |3) the matrix S takes the form (I^ = Sn) 



1 





s 3 / 



(75) 
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where 



Si = cr(s + r - r') + a(r) + cr(r') , 

E 2 = 2cr(s + r - r') + er(s) = C 2 cr(s + r - r') + cr(s) . 



(76) 



Here we show an explicit dependence on the Casimir operator for the large- N c eigenstate 

(77) 



L>f Up 



As a matter of fact, at large N c the quark and gluon colors in the sextet and 15-plet 
states become decorrelated, so that C 6 = C15 = Ca + Cp and 

C A 



C 2 



c f 



(n 



To the leading order in the 1/N C expansion, the state |3) is not excited by single- 
gluon exchange from the initial state |1) = |33). This decoupling is obvious also from 
the projection onto the final states (fTTj) . which at large N c reads 



J2(x\ = £ v dim (^)(^l = v^^ 1 ! + (v^) 3 <2| . 

X R 



(79) 



In the new basis the non-Abelian intranuclear evolution of the 4-body qgqg' state 
becomes the two-channel problem. Expansion over the final states takes the form 



£<X|S[6,E]|1> = ViV c (i|exp 



--ET 
2 



l) + (v/iV c ) 3 (2|exp 



--ET 
2 



(80) 



To the leading order in N c , matrix element of the S-matrix in the first term equals 



(1| exp 



-ET 
2 



|1) = exp 



— EiT 
2 



<^-i[ ff (s) + <7(r) + a(r')]T^ . (81) 



Making use of the Sylvester expansion technique used in for the second matrix 
element one finds 



(2 1 exp 



1 



ET 



|i) = n 



1 exp 



exp 



4e 2 t 



The integral representation of Ref. 



a, 



exp 



exp 



iEoT 



E 2 — Ei 



1 r 1 

-T / d/3 exp 

2 Jo 



E 2 — Ei 



--(/3E 1 + (1-/3)E 2 )T 



(82) 



(83) 
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makes explicit a decomposition into the Initial State and Final State distortions described 
by the cross sections Ei and £2, respectively. Our final result for the sum over final states 
reads 



£<X|S[&,£]|1> = £<X|exp 



x 



X 



-£T 



|1> 



'N c {exp[--[(r(s + r - r') + a(r) + a(r')]T 



, n-T /" 1 rf/3exp[-^(/3£ 1 + (l-/3)£ 2 )T 



+ 



4) 



10-2 

The systematic approach to perturbation 1/N C expansion in has been developed in 
Pi on an example of quark-antiquark dijets in DIS. Its extension to quark-gluon dijets 
is straightforward, we will not dwell into that in this communication. 



V. THE LINEAR /cx-FACTORIZATION FOR DIJETS FORM THE FREE NU- 
CLEON TARGET 



The S-matrices in the master formula (j!4j) depend only on the dipole parameters 
s, r, r'. In the case of the free nucleon target one can integrate over the overall impact 
parameter and represent the integrand of Eq. (|T4*|) in terms of the dipole cross sections: 



2jd 2 bY,(X\ I S§Jb[ } b' g7 b q ,b g )+Sf r (b',b) 



a q*qg + a q*qg' ~ + 

Ca 
C F 



dim(6) 



£21 + 



dim(15) , 



cr(s + r — zr') + <t(s + zr — r' 



\ dim(3)^" 1 ' \\ dim(3) ^ 31 

cr(s + r — r') — cr(s + zr — zr 1 )] 



1 



C A 



[a(s — zr') + er(s + zr) — a(s) — cr(s + zr — zr')} 



+ 77^ 



15) 



Now we apply the fc^-factorizaton representation for the free-nucleon dipole cross section. 
For instance, one readily finds 



f2 = J d 2 H,f(H,)[l — exp(iKr)][l — exp(— ikt')] exp(i«;s) . 



(86) 
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The momentum space wave function of the qg Fock state of the physical quark is 
defined by the Fourier transform 

ty(z,p) = J d 2 r^(z,r)exp(-ipr) . (87) 

We discuss the cross sections averaged over the helicities of the incident parton and 
summed over helicities of the final-state partons. Then ty(z,p) would always enter in 

nn 

combinations of the form [6, 18] 

\*{z,p) - *{z,p- k)\ 2 = 2N c a s ((Q*) 2 )P 9q (z 
where P gq (z) is the familiar splitting function, 



p p — H \ 



p 2 + e 2 (p — n) 2 + e 2 



P gq (z) = C F 1 + {1 z Z) \ (89) 

and, neglecting the mass of the incident light quark, 

e 2 = z(l - z)(Q*) 2 , (90) 

where (Q*) 2 = {p*) 2 is the virtuality of the incident quark q*, given by the square of its 
transverse momentum in the projectile hadron. If e 2 is negligible small compared to p 2 , 
then one can use the large-p approximation, 



\ 2 2 
P p — H \ K, 



p 2 (p — k) 2 J P 2 (p — k) 2 ' 



(91) 



and it is worth to recall the emerging exact factorization of longitudinal and transverse 
momentum dependencies which is a well known feature of the high energy limit. 
Then the master formula for the free-nucleon cross section takes the form 

da N (q* -> qg) _ 1_ 

dzd 2 p q d 2 p g ~ 2(2vr) 

x J d 2 sd 2 rd 2 r' exp[— iAs — ipr + ipr'} exp(iK,s)^l(z, r)ty*(z, r') 



J d 2 Kf{K) 



f C A 

C A 



F 



[1 — exp(zK7*)][l — exp(— kr') 



+ — —[exp(izKr) — exp(?Kr)][exp(— iznr 1 ) — exp(— inr')] 
[1 — exp(izK,r)][l — exp(— iznr')] 1 



iV c 2 -l 
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' r /(A)(^|¥(*,p)-¥(p 



2(27r) 2JV y [2C F 
+ £L\V(z,p- A) - *(p - zA)\ 2 - ^L_|*(^p) - «r(p - zA)| 2 1 



(92) 



A direct comparison shows that the dijet spectrum is precisely the differential form 
of the inclusive single gluon spectrum from the excitation q* — > which was derived 
in (|. The reason emphasized in 16] is that the excitation q* — > qg proceeds via one- 
gluon exchange and the acoplanarity momentum is precisely the transverse momentum 
of the exchanged gluon. Remarkably, the color dipole structure of the integrand of the 
dijet cross section only differs from the one for the single-jet spectrum by the shift of 
arguments of all the dipole cross sections by s. 

The free-nucleon cross-section is a linear functional of the unintegrated gluon density. 
Then, with certain reservations on the region of soft A, the acoplanarity distribution 
is a direct probe of f(x, A). First, on the pQCD side, the unintegrated gluon density 
f(x, A) is well-defined only for sufficiently large momenta A above the soft scale. Second, 
from the practical point of view, any definition of the jet momentum has an intrinsic 
uncertainty with whether the soft hadron belongs to the jet or to the underlying soft 
event, so that experimentally the acoplanarity momentum is well-defined only when it 
is above the soft scale. 



VI. THE NONLINEAR /^-FACTORIZATION FOR THE DIJET PRODUC- 
TION OFF NUCLEI 

A. The color-dipole representation at large N c 

The final Fourier representation for the leading term of the large-iV c expansion for 
the dijet cross section per unit area in the impact parameter space reads 

da(q* -> qg) 1 



d 2 sd 2 rd 2 r' 

d 2 bdzd 2 Ad 2 p (2vr) 4 ' 

x exp[— 2 As — ipr + ipr'}^(z, r)^*(z, r') 



~n-T(b) [ t^exp 
2 Jo 



--[/?£! + (1-/?)£ 2 ]T(6) 
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+ exp 
+ exp 

— exp 

— exp 



zr + zr)T{b) 



-~[<r(s + r-r') + a(r) + a(r')]T(b) 
1 f 

--[a(r)+a(s + r-zr'))T(b) 
--[o-(r') + a(s -r' + zr)]T(b) 



(93) 



Recall that the first term, oc Q, describes the excitation from the color-triplet dipole to 
sextet and 15-plet dipole states. Note, how the large- N c suppression of the off-diagonal 
matrix element £12 is compensated for by a large number of final states in the higher 
representations. At large N c , once the sextet and 15-plet states have been excited, 
their de-excitation back to the triplet state is suppressed and the further intranuclear 
evolution consists of the color rotations within the higher representations. The remaining 
four terms in describe the rotations within the color triplet states. 



B. Unintegrated collective nuclear glue and isolation of initial state distortions 

The transformation from the color-dipole to the momentum representation is fur- 
nished by the /cx-factorization formula ()19|) and its generalization to the nuclear target. 
For the latter we adopt the collective nuclear unintegrated gluon density per unit area in 
he impact parameter plane, <p(b, x, k), as defined in terms of the nuclear profile function 



4. 



3CJ: 



r[6, a(x, r)] = 1 — exp 



~a(x,r)T{b) 



d 2 K(f)(b, x, k) 



(94) 



The utility of <p(b, x, k) stems from the observation that the driving term of small-x 
nuclear structure functions, the amplitude of coherent diffractive production of dijets 
off nuclei and the single-quark spectrum from the 7* — > qq excitation off a nucleus all 
take the familiar fc^-factorization form in terms of <j)(b,x,K,). The so defined collective 
nuclear glue 0(b, ) satisfies the sum rule 

d 2 K(fi(b, x, k) = 1 - S[6, a (x)] , (95) 

where (Jq(x) is the dipole cross section for large dipoles. The multiple-scattering ex- 
pansion of 0(6, x, k) in terms of the collective glue of j-overlapping nucleons in the 
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Lorentz-contracted nucleus and its nuclear shadowing and antishadowing properties are 
found inflflQHand need not be repeated here. We on, cite the form nla for the 



so-called saturation scale 



4tt 2 

Ql(b,x) w — a s (Q 2 A )G(x,Q 2 A )T(b) 



(96) 



and reiterate that at a large saturation scale <p(b, x, k) is well-defined not only for per- 
turbative values of n 2 below Q\(b,x), its continuation to the soft region is also stable. 
To this end we recall that although <Jq(x) enters the multiple-scattering expansion for 
<p(b,x, k), the final form of <j)(b,x, k) is exclusively controlled by Q\(b,x) and does not 
depend on the auxiliary soft parameter <7o(x) |4]. We also note, that the nuclear pro- 
file function satisfies the s-channel unitarity bound for the partial waves of the dipole- 
nucleus scattering, T[b,cr(x,r)] < 1, while the partial wave of the impulse approxi- 
mation (IA) overshoots the s-channel unitarity bound for sufficiently heavy nucleus, 
r( /A ) [b, a(x, r)] = \g{x, r)T{b) > 1. As such, the unintegrated collective nuclear gluon 
density <p(b, x, k) defined by Eq. ()94)1 unitarizes the density of partons in a Lorentz- 
contracted ultrarelativistic nucleus. 
Still another convenient quantity is 



$(6, x, k) = S[b, o- (x)}5(k) + (j)(b, x, k) 



(97) 



in terms of which 



exp 



~-a(x,r)T(b) 



J d 2 K&(b, x, k) exp(«Kr) 



We shall also encounter the collective glue for a slice < (3 < 1 of the nucleus: 



cxp 



-~Pa(x,r)T(b) 



d 2 K&((3; 6, x, k) exp(inr) 



(9h 



(99) 



and the intranuclear attenuation-distorted wave function in the dipole and momentum 
representations, 

1 



ty((3,x;z,r) = ^(z, r) exp 



-3a(x,r)T(b) 



*(/3, X] z } p) = J d 2 r^(/3; z, r) exp(-ipr) = J d 2 K^(z, p - «)$(/3; 6, x, k) . (100) 

Hereafter, unless it may cause a confusion, we suppress the variable x in gluon densities, 
dipole cross sections and distorted wave functions. 
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C. Excitation of color-triplet quark-gluon dipoles 

First, we rewrite the last four terms in the integrand of ()93j) in terms of the distorted 
wave functions. Then we make use of the Fourier representation (|98|). (jlOOj) : 

da{q* -> qg) 



d 2 bdzd 2 Ad 2 p 



— - [ d 2 sd 2 rd 2 r' exp[— 2 As — ipr + ipr'] 

7T) 4 J 



(2tt 

| ^(l;z,r)^*(l;z,r')exp 
+ #(2, r)#*(z, r')exp 



-tr(s + r - r')T(b) 



- ^(l;z,r)^*(z,r')exp 



-<t(s - zr' + zr)T(b) 
—a(s + r - zr')T(b) 

- *(z, r)**(l;z,r')exp -^cr(s - r' + zr)T(b) 

= \ / d 2 sd 2 rd 2 r'd 2 K$(b, k) exp[-iAs - ipr + ipr'] 

(z7r) 4 J 

l *(1; z, r)\P*(l; 2, r') exp[iK(s + r - r')] 
+ \I/(z, r)^*(z, r') exp[iK,(s — zr' + zr)\ 

- ^(1; z, r)**(z, r') exp[iK;(s + r - zr')] 

- #(z, r)**(l; z, r') exp[iAt(s - r' + zr)] | 

1 



(2tt) 2 
1 

1 

(2^ 



$(&, A) |*(l;z,p- A) - #(z,p-zA)| 
0(6, A) |tt(l; z,p — A) — #(z,p - zA)| 
6(A)\^(l;z,p)-^(z,p)\ 2 S[b,a (x)] 



(101) 



Now recall p} that the amplitude of the coherent diffractive excitation qA — > (g^)^4 is 
precisely proportional to 



*(z,p) - #(l;z,p) = J d 2 r^(z,r) |l - exp 



~a(r)T(6) 



exp [—zpr J 



(102) 



so that the last term in (jlOlj) describes the coherent diffractive production of dijets. In 
the approximation of very large nucleus the diffractive dijets are produced exactly back- 
to-back. For finite nuclei instead of the delta-function 5(A) one finds the sharp peak of 
the width A 2 < 1/R\ which is described by the form factor of the nucleus, the details 
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are found in |8[ and must not be repeated here. The former term describes inelastic, 
incoherent production of color-triplet qg states. 



D. The contribution from sextet and 15-plet final states 



The evaluation of the contribution from the excitation of higher color representations 
in (j^Hj) proceeds as follows. First, we make use of the integral representation (jSfij) for 
the off-diagonal cross section. Second, keeping an explicit dependence on the Casimir 
operators 6*6,15, we have 







dp exp 



10 

x exp 

= / 

Jo 

X 



[1 - 0)E 2 )T(6) 



CM V 



[l-P)[C 2 a(s 



--Pa(r)T(b) 



a(r')]T(b) 
a(s)]T(b) 



exp 



l -Pa{r')T{b) 



dp exp 

f 1 

2 

dp exp 

J d 2 K$(P; b, k 3 ) exp[zK 3 (s + r - r')) 
x J d 2 K 2 $(C 2 (l - P); b, k 2 ) exp[iK 2 (s + r- r')} 
x J d 2 Ki$(l — P; b, Ki) expfzKis] 



(103) 



In this decomposition we keep the dipole form of the two attenuation factors S[b, P<j(r)] 
and S[b, P<j(r')]. They describe the coherent intranuclear distortion of the color-triplet 
quark-gluon dipole before the excitation into the sextet and 15-plet representations at 
the depth P from the front face of the nucleus. The way to handle these distortion 
factors has already been clarified above. Note, that in contrast to the quark- ant iquark 
dijet production in DIS off nuclei, both the ISI and FSI distortion factors depend on the 
dipole parameter s and explicitly contribute to the acoplanarity distribution. 

Combining together (JHtjj). (jlOOJ) and (jl(J3|) we obtain the dijet spectrum from the 
excitation of the sextet and 15-plet dipoles 



da(q* -> qg) 



6+15 



(2n 



r2 T {b) j o d p 



d 2 bdzd 2 Ad 2 p 
x J d 2 Kd 2 Kid 2 K 2 d 2 K35(K + K\ + n 2 + n 3 — A) 
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x - (3; b, Kl )$(C 2 (l - 0); b, k 2 )$(/3; 6, /c 3 ) 

x |*(/3;«,p-«2-K 3 ) -^(/?;2,p-k 2 -k 3 -k)| 2 (104) 

The acoplanarity momentum A manifestly receives four distinct contributions which can 
be classified as follows. The excitation from the color-triplet to the sextet and 15-plet 
states by single-gluon exchange with one of the nucleons of the nucleus contributes the 
transverse momentum k. The momentum k 3 comes from the ISI of the incident quark, 
the FSI of the qg dipole in the sextet and 15-plet representations contributes Kx and k 2 . 

The emergence of the collective nuclear glue $(6*2(1 — b, k 2 ) in the integrand of 
(jl()4j) is not accidental. While (1 — (3) is a thickness of the slice of the nuclear matter 
traversed by the sextet and 15-plet qg dipoles, the factor C 2 derives from the Casimir 
operators of higher representations, see Eq. (|77|). That is one more illustration of our 
point jj, 0| that the collective gluon field of the nucleus cannot be described by a single 
density function, it is a density matrix in the space of color representations. In the 
considered large- N c approximation, C 2 = Ca/Cf and $((1 — 0)C a/ Cp\ b, n 2 ) is precisely 
the collective nuclear glue defined in terms of the color-singlet gluon-gluon dipole. 

The ISI and FSI distortions can partly be combined taking the convolution ^ 

J d 2 K 3 d 2 K 2 $(C 2 (l -/3);b, «2)$03; 6, k 3 )8(k -k 2 - k 3 ) = $(/? + C 2 (l - (3); b, k) . 

(105) 

which is also obvious from the color-dipole form in ()103j) . 



E. Nonlinear /^-factorization for dijets: the universality classes 

1. Quark-gluon vs. quark- antiquark dijets 



After application of the convolution (JTUoJ), the final result for the inclusive dijet spec- 
trum takes the form 

(2n) 2 da A (q* ^ qg) 1, * 



-T{b) J d{3 J d 2 K 1 d 2 Kf{x, k) 



d 2 bdzd 2 pd 2 A 
x$(l - p, b, x, A - k x - k)$(J3 + C 2 (l -P),b, x, k x ) 

2 

x <Sf(p;z,p- Kx) - V(p;z,p- Kx - k) 
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+0(6, x, A) tf(l; z, p - A) - V(z, p - zA) 
+5(A)S[b,a Q {x)) y(l ]Z ,p)-V(z,p) 2 . 



(106) 



which must be compared to the large- N c version of the free-nucleon cross section (|92|) . 

The free-nucleon cross-section is a linear functional of the unintegrated gluon density. 
The fcj_-factorization properties of the nuclear cross section are much more complicated. 
At this point, it is instructive to discuss ()106|) in conjunction with the quark- ant iquark 
dijet spectrum in DIS [4] and gluon-nucleus collisions [l(|. The spectrum of dijets in 
DIS equals 



where the first term describes the excitation of the color-dipole from the lower (color- 
singlet) to higher (octet) representation, whereas the second term is the contribution 
from coherent diffractive excitation. The spectrum of the quark-antiquark dijets in gA 
collisions is of the form 



Now we can identify the four universality classes of the nonlinear fc^-factorization which 
differ by the pattern of transitions between the initial and final state color multiplets. 
They describe the leading transitions in the large-iV c approximation, the higher order 
excitation and regenerationprocesses result in still higher nonlinearity in gluon densities, 
the examples are found in 4|. 

2. Excitation of higher color representations from partons in the lower representations 

Excitation of color-octet states in DIS, and of sextet and 15-plet states in qA inter- 
actions, belong to this universality class. The two reactions have much similarity. In 




x/(k)$(1 - P, b, A - k x - - P, b, Kl ) 

2 

x ty(p',z,p- «i) - V(P;z,p- k x - k) 
+5{A) *(l;z,p) -*(z,p) 2 . 



(107) 




x\^(z,p_ - k) - ^(z,p_ - zA)\ 2 . 



(108) 
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both cases the nonlinear fcj_-factorization formulas contain the free-nucleon gluon density 
f(x, n), which describes the transition from the qg color dipole from the lower - triplet 
for qg and singlet for DIS - to higher - sextet and 15-plet for qg and octet in DIS - color 
dipoles. In both cases, the number of states in higher representations is by the factor 
N% larger than in the lower representation. In qq excitation in DIS the corresponding 
contribution to the dijet spectrum is the fifth order functional of gluon densities. In 
the qg case it is the sixth order functional of gluon densities, only after the application 
of the convolution (jl()5j) it takes the form of the fifth order functional. Two powers of 
the collective nuclear glue enter implicitly via the coherent ISI distortions of the wave 
function \[f(j3]z,p) in the slice of the nuclear matter before excitation of color dipoles 
in the higher representation, two more powers of the collective nuclear glue describe the 
ISI and FSI broadening of the acoplanarity distribution. 

The principal difference between DIS and qA interactions is in the nuclear thickness 
dependence of the distortion factors. Namely, the factor 

$((1 — P), b,A — Ki — k)$((1 -/3),b, Kl ) 

in DIS is the symmetric function of the nuclear gluon momenta k,i and k 2 = A — K\ — k, 
which flow from the nucleus to the quark and antiquark (or vice versa), respectively. It 
describes equal, and uncorrelated, distortion of the outgoing quark and antiquark waves 
by pure FSI. The independence of the two distortion factors is a feature of the large N c 
approximation. For qg dijets in qA collisions the distortion factor 

$(1 - (3, b, « 2 )$(C 2 (1 ~P)+P, b, ki) 

is an asymmetric one. The first source of the asymmetry is the non-singlet color charge 
of the projectile parton. The second source is that the two partons in the final state 
belong to different color representations. This is best seen from in the overall distortion 
factor in (fTTH|) . 

$(# b, «3)$(C 2 (1 -P);b, k 2 )$(1 - (3- b, Kl ) , 

before taking the convolution (|1U5|) . The FSI distortions in the slice (1 — f3) of the 
nucleus are given by the two last factors, of which $(1 — /3; b, K\) is a broadening due 
to final-state rescatterings of the quark. Because C 2 = Ca/Cf, see Eq. (JZ8J), the second 
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FSI factor, $(C 2 (1 — /?); 6, k 2 ), corresponds to the FSI distortion of exactly the outgoing 
gluon wave. To the large- N c approximation the rescatterings of the quark and gluon are 
uncorrelated. 

The coherent ISI distortion of the wave functions in DIS and qA collisions is identical. 
However, in qA collisions this coherent distortion is accompanied by an incoherent ISI 
distortions of the incident quark wave described by b, K3). In DIS the incoherent 
ISI distortions are absent because the photon is a color-singlet particle. We can anti- 
cipate that gluon-nucleus collisions with excitation of gluon-gluon dijets in higher color 
representations will belong to this universality class. 

3. Excitation of final state dipoles in exactly the same color state as the incident parton: 
coherent diffraction 

To this universality class belong the exactly back-to-back dijets. Another experimen- 
tal signature of the coherent diffraction is a retention of the target nucleus in the ground 
state and large rapidity gap between the hadronic debris of the diffractive dijet and the 
recoil nucleus. It is most important for DIS where coherent diffraction dissociation of 

n 

the photon into qq dijets makes for heavy nuclei ~ 50% of the total DIS rate |3J|. The 
origin of the coherent diffraction is a coherent nuclear distortion of the wave function of 
the qq Fock state over the whole thickness of the nucleus. 

In the coherent diffractive excitation of qg dipoles in qA collisions the qg dipole must 
propagate in exactly the same color state as the incident quark. The nuclear suppression 
factor S[b, ao{x)} has the meaning of 

S[b, f 7 (x)] = (s[6,i ( To(x)]) 2 (109) 

and the factor S[b, |o"o(a;)] in the diffractive amplitude corresponds to the intranuclear 
attenuation of the quark wave with the total cross section 

v qN = 2°"o(^) • (110) 

Coherent diffractive excitation of color-octet gluon-gluon dijets in gluon-nucleus collisions 
is expected to exhibit similar properties. 
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Coherent diffract ive excitation of QQ dipoles in gA collisions is allowed, but it is 
suppressed at large N c by the condition that the QQ dipole must propagate in exactly 
the same color state as the incident gluon. 

4- Incoherent excitation of final state dipoles in the same lower color representation as the 
incident parton 

An example of this universality class is an inelastic excitation of color-triplet qg states 
in qA collisions followed by a color excitation of the target. Here both the incident 
parton and dijet belong to the fundamental, i.e., lower, representation of SU(N C ). The 
intranuclear evolution of such a dipole is confined to rotations within the color-triplet 
state. This contribution is not suppressed at large N c . The dijet cross section for this 
universality class looks like satisfying the linear /cj_-factorization in terms of <p(b, x, A). 
But this is not the case: one of the wave functions, z,p g ), is coherently distorted 
over the whole thickness of the nucleus, so that this contribution is a cubic functional of 
the collective nuclear glue. 

We can anticipate that gluon-nucleus collisions with excitation of color-octet gluon- 
gluon dijets will belong to this universality class, although one has to account for the 
existence of the two, F-coupled and D-coupled, octet states. 

Although superficially it looks like a subclass of this universality class, the coherent 
diffraction is a distinct class for the property of the exact back-to-back dijets and the 
rapidity gap between the dijet and the recoil nucleus in the ground state. 

5. Excitation of final state dipoles in the same higher color representation as the incident 
parton 

In the realm of QCD with gluons in the adjoint representation and quarks in the 
fundamental representation, this universality class consists of the quark- ant iquark di- 
jets in gluon-nucleus collisions. Only in this case the initial parton (gluon) belongs to 
the higher (octet) color multiplet of the final QQ state. At large N c , the intranuclear 
evolution of QQ will consist of color rotations within the space of color-octet states. 



39 



The de-excitation from the color-octet to color-singlet QQ dipoles is suppressed at large 
N c . Consequently, the non-Abelian evolution of the QQQ'Q 1 state becomes the single 
channel problem. The coherent diffraction excitation, in which the initial and final color 
states must be identical, is likewise suppressed. The emerging pattern of quadratic non- 
linearity can be related to the large- N c gluon behaving like the color-uncorrelated quark 
and antiquark. 

The above classification exhausts reactions caused by incident photons, quarks and 
gluons. However, technically all the universality classes have a much broader basis. 
Indeed, instead of an incident gluon one can think of the projectile which is a compact 
lump of many partons in the highest possible color representation. For instance, compact 
diquarks in the proton can be viewed as sextet partons. 



6. Is an experimental separation of events belonging to different universality classes possi- 
ble? 

We reiterate that for all the universality classes their separate contributions to the 
dijet cross section are infrared-safe quantities. Coherent diffraction has distinct sig- 
natures and the experimental separation of events from this universality class is not 
a problem. Production of very forward dijets in proton-nucleus collisions evidently 
tags quark-nucleus collisions. Production of open charm in the proton hemisphere of 
of proton-nucleus collisions tags gluon-nucleus collisions. Incoherent processes belonging 
to different universality classes are characterized by distinct color charge of the hard dijet 
and this distinction is well defined at the parton level. Translating the cross-talk be- 
tween color charges in the dijet, the spectator partons of the proton and the color-excited 
nucleus remnant into properties of hadronic final states can only be done within nonper- 
turbative hadronization models. As an example we cite the impact of color reconnection 
effects on the flow of slow hadrons and the accuracy of the mass determination in 
e + e~ annihilation (jsjj, for the review see js^). 
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F. The impulse approximation 

In the impulse approximation (IA) one only has to keep the terms linear in T(b). 
The transition to the IA is best seen in the color-dipole representation (|93p. Recall, that 
our formulas for nuclear cross section were derived in the large- N c approximation. Here 
the first term, the contribution from the sextet and 15-plet final states, is already linear 
in T(b) and one must put the attenuation factors equal to unity. The remaining four 
exponentials must be expanded to terms linear in T(b). Then one would find precisely 
the large- N c version of Eq. (|85j) times T(b). The integration over impact parameters 
gives / d 2 bT(b) = A. Such a comparison does not expose the role of coherent diffraction 
and we revisit the issue in the momentum representation. 

We start with the sextet and 15-plet contribution in ()106|) . It already contains 
the factor T(b). Consequently, one must neglect ISI distortions in the wave function, 
y(P;p) =>^(p), and take 

$(1 - (3, b,x, A - «! - k)$(2 - 0, b, x, ki) = 8{A -k x - k)5{k x ) . (Ill) 

This way one would recover the first term in the rhs of Eq. (|92jl . In the contribution 
from the excitation of the triplet dipoles one must neglect the distortion of the wave 
function and take 

4>(b,X,K) = ^T(b)f(x,K,). (112) 

The second term in the rhs of Eq. (}9*2~|) is recovered. Finally, according to Eq. (jlU2|) the 
diffractive amplitude starts with the term linear in T(b). Consequently, the coherent 
diffractive contribution to the dijet cross spectrum starts with the terms oc T 2 {b) and 
vanishes in the IA. 

VII. NUCLEAR BROADENING OF THE ACOPLANARITY DISTRIBUTION 

The nuclear broadening of the acoplanarity distribution of hard quark-gluon dijets 
from qA collisions is somewhat different from the broadening of quark-antiquark jets in 
DIS and now we comment on those differences. 
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A. Coherent diffractive contribution 



The first striking difference is in the role of the coherent diffractive production. It gives 
exactly back-to-back dijets. In the considered approximation of single-gluon exchange in 
the t-channel diffractive production off the free-nucleon target vanishes. Experimentally, 



3- 
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at HERA energies a fraction of DIS which is diffractive does not exceed 10% 
contrast to that, in DIS off heavy nuclei a fraction of coherent diffraction was shown 
to be as large as ~ 50% [sfj]. The existence of coherent diffractive mechanism in the 
quark-nucleus collisions is interesting by itself. From the practical point of view, it is 
suppressed by nuclear absorption and is marginal. 

B. Excitation of the color-triplet states 

Inelastic excitation of color-triplet dipoles is a specific feature of qA collisions in the 
sense that it has no counterpart in DIS. One must compare 

<f)(b,x, A)|tf (1; z,p - A) - V{z,p - zA)\ 2 (113) 

with its IA form 

~T(b)f(A)\*(z,p- A) - ¥(p - zA)\ 2 . (114) 

The first striking distinction is that that for the free-nucleon target the contribution 
of this process vanishes at z — > 1, when the incident quark's momentum is transferred 
entirely to the forward gluon jet. For the nuclear target this is not the case because one of 
the wave functions in (jllHJ) is the nuclear-distorted one. Because p — A = — p q , it takes 

2 

the form 4>(b,x, A) ^(l; z, p ) — ^/(z,p q ) ; as a function of the quark-jet momentum, it 
is reminiscent of the coherent diffractive contribution, but the acoplanarity momentum 
distribution is given by the unintegrated nuclear gluon density <p(b,x, A). Hereafter we 
consider the case of finite (1 — z). 

A comprehensive discussion of nuclear properties of the ratio 
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is found in |4|, It is nuclear-shadowed, R g (b,A) < 1, for A < Q\{b) and it exhibits 
antishadowing property, R g (b, A) > 1 in a broad region of A 2 > Q\(b) . The maximum 
value of R g {b, A) is reached at a value of A 2 which is larger than Q\(b) by a large 
numerical factor. 

Now we turn to distortions of the wave function. We are interested in hard dijets. 
If the incident quark is a valence quark of the proton, its transverse momentum and 
virtuality have the hadronic scale and can be neglected. For hard jets 

*(z,p)cx4 (H6) 
and, upon averaging over the azimuthal angle <p of the gluon momentum n, 

(*(z,p-K))„<x^e(p 2 - K 2 ). (117) 
Consequently, the wave function distortion factor equals 



P*(b,Z,p)= ; ' = / d*K$(b,K) = l- / d 2 K^(b^). (118) 

ty{z,p) J Jp 2 



For the weakly virtual incident quark it does not depend on z. For hard jets, p 2 > Q 2 (b), 
the remaining integral ()118j) can be evaluated following the analysis of the Cronin effect 
in 0. Namely, here we can use the leading-twist approximation, 

$(b, K ) = ir(6)/( K ), (ii9) 

and the definition (|2Uj) with the result 

<!L = 1 -pJb,z,p) = / d 2 K$(b, x, k) 

Jp 2 

^ 2n 2 T(b)a s (p 2 ) dG(x,p 2 ) = I Q|(6) a s (p 2 ) dG(x,p 2 ) 

N cP 2 dlogp 2 2 p 2 a s {Q\)G{x,Q 2 A ) d\ogp 2 ' l ' 

It is important that 5^ is a manifestly positive valued quantity. It has a form similar to, 
but is numerically smaller than, the nuclear higher twist correction to <p(b, x, k). 

In Fig. ED we show the numerical results for the wave-function distortion factor for 
the gold nucleus at several values of the optical thickness u{b) = ~a (x)T(b). At this 
point one needs to pay a due attention to an explicit dependence on the QCD running 
coupling Qtg(r) on the small dipole size r in Eq. (|2()jl. The discussion of its impact is 
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FIG. 12: The left panel shows the impact-parameter dependence of the optical thickness of the 
gold nucleus for several values of the gluon-jet momentum p. The momentum dependence of 
the wave- function distortion factor p^(b, z,p) for several values of the optical thickness of the 
nucleus is presented in the right panel. 



found in |4j, |£| , in the evaluation of the momentum spectra this running coupling must be 
taken at the largest relevant hard parameter, which in our case is p 2 . Correspondingly, 
in all the formulas for the dijet spectra, the dipole cross section for large dipoles, cr (x), 
must be understood as 

4tt 2 rdK 2 



l\ r J K 



'121 



For this reason, the optical thickness of the nucleus v(b) as a function of the impact 
parameter b, shown in the left panel of Fig. depends on the hard scale - the jet 
momentum. The wave-function distortion factor p^(b,z,p) is shown in the right panel 
of Fig. El The hard regime ()120)) for 5^ sets in at the momenta p > 1 GeV. We reiterate 
that the saturated cross section <Jq(x, 00) is only an auxiliary parameter which does not 
enter directly the observable cross sections - the latter only depend on the saturation 
scale Q\(b), the discussion is found in Ref. ^J. 
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In terms of the distortion factor p^(b,p) one readily finds 
*(l;z,p-A)-*(z,p-z&) 



\^(z,p - A) - $(p-zA)\ 2 
[(1 - z)A - 5+(p - zA)f _ [(1 - z)A + 64 Pq - (1 - z)A)f 



122) 



(1-^) 2 A 2 (1-^) 2 A 2 

The overall nuclear modification factor, the ratio of the nuclear, (J113I) . and free-nucleon, 
(jll4j) . target contributions, is a product 

Rf /N (b,p,A) = R g (b, A)R^(b,p, A) (123) 

Here R g (b, A) does not depend on the jet momentum p except for the weak depen- 
dence through as{p 2 )- Evidently, R^(b, p, A) is azimuthally asymmetric and favors 
A anticollinear to the gluon momentum and collinear to the quark momentum: in the 
back-to-back configuration, the gluon jet tends to have the transverse momentum smaller 
than the quark jet. The dominant contribution to the nuclear dijet cross section comes 
from A 2 ~ Q\(b), and for hard dijets the asymmetry will be weak, of the order of 
~ Qa(P)/p. 

Alternatively, if one keeps the quark transverse momentum fixed and increases the 
target mass number A, i.e., Q\(b) and 5$ thereof, the transverse momentum of the 
away gluon jet will decrease with A. The form of the q — > qg splitting function favors 
production of the gluon jet at rapidities smaller than the quark jet. Then, the above 
correlation between the acoplanarity and quark momenta shall exhibit itself as a nuclear 
suppression of the away jet produced at rapidites smaller than the rapidity of the forward 
trigger jet. The numerical studies of this effect will be reported elsewhere. 

C. Excitation of the sextet and 15-plet jets states 

Here one must compare the contribution to the nuclear dijet spectrum (jlU4j) with its 
IA counterpart 



T(b)da N (p,A) 



dzd 2 pd 2 A 



6+15 



2(2tt 



' T(b)f(A)\^(z,p)-^(p-A)\ 2 . (124) 
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Note, that the nuclear cross section can be cast in the from reminiscent of a triple 
convolution 

da A (q* -> qg) 



6+15 



cPbdzdpdA 

T(b) J d(3 J d 2 Kd 2 Kid 2 K 2 d 2 K?,5{K + Ki + k 2 + K 3 - A) 

da N (p - k 2 - k 3 , k) 



x $(1 - (3; b, Kl )$(C 2 (l -/3);b, k 2 )$(/3; 6, k 3 ) 



• (125) 



6+15 



dzd 2 pd 2 n 

which suggests that at a fixed gluon-jet momentum p, it will be a broader distribution of 
A than the free-nucleon cross section (for the related discussion see jj]). This broadening 
is best seen for hard dijets, p 2 3> A 2 ,Q\(b). Because the dominant contribution comes 
from k, 2 < Q\(b), one can neglect ^2,3 compared to p in the free-nucleon cross section 
in the integrand of (jl25j) . Then the nuclear cross section takes the manifest convolution 
form 



d(T A (q* -> qg) 



6+15 



d 2 bdzdpdA 



(126) 
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The saturation scale for the distribution $(1 + C 2 (l — /3); b, A — k) equals 



[1 + C 2 (l - P)}Q 2 A (b) 



(127) 



and the broadening of the acoplanarity distribution for the quark-gluon dijets is sub- 
stantially stronger than that for the quark- ant iquark dijets in DIS discussed in j^J. 



VIII. THE MONOJETS FROM DIJETS: FRAGMENTATION VS. GENUINE 
DIJETS 

A. Monojets from dijets in the free-nucleon reactions 

In the above discussion we implicitly assumed that the quark and gluon hard jets are 
separated by a large azimuthal angle and the acoplanarity momentum is small compared 
to the jet momenta, A 2 < p 2 , (p — A) 2 . The interesting new situation is encountered 
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when the quark and gluon jets start merging. Specifically, the wave function ^(z, p—zA) 
has a pole when p — zA = 0, i.e., when the gluon and quark are collinear, 



P g = zA, Pq = z q A = A-p= {1- z)A = z q A. (128 

In the vicinity of the pole the qg production cross section has the factorized form 
da N (q* -> qg) 

monojet 



dzd 2 pd 2 A 



-^f(A)\^(z,p-zA)\\ (129) 



Now recall that ^(z,p — zA) is precisely a probability amplitude to find the gluon 
with the momentum = p — zA transversal with respect to the axis of the quark jet 
with the momentum A, and \^(z,p — zA)\ 2 of Eq. (|KK|) is proportional to the familiar 
splitting function P gq (z), which is precisely the driving term of the quark-jet fragmenta- 
tion function. Consequently, the contribution ()129|) must be treated as a fragmentation 
of the scattered quark into the quark and gluon, q' — * qg. The quark pole contribution 
will dominate if 

fci«(p-A) 2 = pJ. (130) 

From the experimental point of view, the corresponding final state is a monojet of the 
transverse momentum A. The transverse momentum of such a monojet will be compen- 
sated by an away jet produced at midrapidity or the nucleus hemisphere of pA collisions. 

In terms of Feynman diagrams of Fig. |21 - for the free-nucleon target one takes 
the single-gluon exchange, - the monojet production is a property of the diagram (c). 
Indeed, the cross section (I129j) is proportional to precisely the differential cross section 
of quasielastic scattering of the projectile quark off the nucleon target - the latter is 
evidently proportional to the unintegrated gluon density of the target proton /(A). The 
two classes of Feynman diagrams in Fig. El (b) and (c), are integral parts of the gauge- 
invariant description of the QCD Bremsstrahlung excitation of the qg state. Still, the 
isolation of the pole contribution from the gauge-invariant combinations 



ty(l;z,p- A) - ty(z,p- zA) = ty(l;z,p g ) —^f(z,p — zA 



2 



2 



in (|92|). and of the monojet contribution to the generic dijet cross section wouldn't conflict 
gauge invariance. In order to conform to the jet-finding algorithms, the production of 
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the quark and gluon within the jet-defining cone must be treated as a fragmentation 
of the monojet; if the azimuthal angle between the quark and gluon is larger than the 
jet-defining angle, the two jets must be viewed as independent ones. The combination 
of the wave functions, which enters the excitation of the sextet and 15-plet final states, 
see Eq. (II 24(1 . has the form 



and is finite for all orientations of the quark and gluon jets. 

The quark-tagged pQCD gluon Bremsstrahlung considered here is already the higher 
order process, the lowest order pQCD process in qN interaction is the radiationless 
quasielastic scattering of the quark. Naive application of fragmentation q' — > qg to this 
lowest order process would evidently lead to a double counting, because the fragmen- 
tation is manifestly a monojet part of our dijet cross section. The integration over the 
gluon momentum k± in the inclusive cross section would yield the familiar collinear log- 
arithm, which must be reabsorbed into the definition of the fragmentation function at 
the starting scale. Simultaneously, one must include the virtual radiative correction to 
the radiationless quasielastic scattering of the incident quark off the target nucleon. The 
treatment of these virtual corrections to quasielastic scattering and elimination of double 
counting go beyond the scope of the present study and will be addressed elsewhere. We 
only want to comment that if one would insist on the description of monojets in terms of 
the fragmentation of the quark, then the interplay of the virtual correction to the radia- 
tionless quasielastic scattering and of the collinear logarithm in the monojet component 
of the the dijet cross section may entail a departure of the fragmentation function from 
that defined in the e + e— annihilation. 

B. Monojets from dijets off a nuclear target 

The presence of the monojet pole ()128|) in the nuclear dijet cross section ()106|) is 
manifest: 



ty{z,p) - ^{z,p- A)| 2 oc 




daAf -> qg) 

d 2 bdzdpdA 



monojet 



2(2vr) 



1 



-T(b)<f>{b,x,A)*{z,p-zA) . (131) 
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It factorizes precisely as the free-nucleon cross section: the differential cross section of 
quasielastic quark-nucleus scattering, proportional to the unintegrated collective gluon 
density of the nucleus, times the fragmentation of the scattered quark to the gluon and 
quark given by \^(z,p — zA)\ 2 , which does not depend on the target. However, the vir- 
tual radiative correction to the radiationless quasielastic scattering of the incident quark 
off the target nucleus and the elimination of double counting are likely to depend on 
the acoplanarity momentum A and the shape of the collective nuclear glue (p(b, x, A). 
Should this be the case, such a dependence could be reinterpreted as a nuclear modifi- 
cation of the fragmentation function; this issue will be addressed elsewhere. 

As it was the case for the free-nucleon target, excitation of the sextet and 15-plet 
final states is free of the monojet singularities. To be more precise, the wave-function 
singularities in the integrand of the sextet and 15-plet contribution to ()106)) occur in the 
intermediate state, at p — k,± — n = and p — k± = 0. However, they are integrated 
out in the observed dijet cross section. It is still instructive to look at the effect of these 
singularities in the monojet kinematics A 2 ^> p 2 >,Q\(b). 

Consider first the contribution from the intermediate pole at p — K\ = 0. The relevant 
Ki integrations are of the form 

J d 2 K 1 d 2 Kf(x, - P, b, x, A - «i - k)$(P + C 2 (l - P), b, x, «i) 

x |tf(/3;*,p-/ci)| 2 

= ${/3 + C 2 {l-P);b,x,p) J P d 2 k\^(f3;z,k)\ 2 

x J d 2 Kf(x,K)$(l- {3,b,x,A- p- k) (132) 

For the considered hard jets 

$(1 - 0, b, x, A - p - k) = i(l - P)T(b)f(A -p-k) (133) 

and the convolution in ()132|) equals jj, |^ 

J d 2 Kf(x, - /3, b, x, A - p - k) = (1 - P)T(b)f{A - p) . (134) 

The resulting contribution from the intermediate pole of the wave function at p — Hi = 
is proportional to 

T 2 (b)f(A - p)f{p)P 9q {z) = T 2 (b)f(p g )f(p q )P gq (z) (135) 
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and has the form of the product of the differential cross sections of independent quasielas- 
tic scattering of the quark and gluon fragments of the incident quark. It does not depend 
on the azimuthal angle between the quark and gluon jets at all, and has no collinear 
singularity. A similar situation has been found to occur in our previous study of the 
production of hard quark-antiquark dijets in ttA collisions |l3(. The contribution from 
the pole at p — K\ — k, = is entirely similar. 

IX. CONCLUSIONS 

We presented a derivation of nuclear modifications of the quark-gluon production in 
quark-nucleus collisions. Our principal result is the nonlinear fc^-factorization relation 
(jlOfij) . The derived dijet cross section can be decomposed into three major contributions. 
The excitation of qg dijets in higher - sextet and 15-plet - color representations gives rise 
to the sixth order nonlinearity in gluon fields, compared to the fifth order nonlinearity 
for qq dijets in DIS. A part of the nonlinearity comes from the free-nucleon gluon density 
which emerges in all instances of excitation of higher color representations (see also 
the related discussion of the 1/(N% — 1) expansion in Ref. Q]). The matrix elements of 
transitions from lower to higher color representations are suppressed at large N c , but this 
suppression is compensated for by the large number of states in higher representations. 
The coherent diffraction, in which the final dipole is produced in exactly the same color 
state as the incident quark, is not suppressed by large N c , but because of the color- 
nonsinglet incident partons the diffractive contribution is suppressed by an overall nuclear 
attenuation and will only come from collisions at the diffuse edge of a nucleus. A new 
feature of qA collisions in contrast to DIS is inelastic production of qg states in the 
same color representation as the incident parton. Such color rotations within the same 
representation are not suppressed at large N c . This contribution has the form which 
superficially looks like satisfying the linear /c^-factorization in terms of the collective 
nuclear gluon density. However, it contains the nuclear-distorted wave function of the 
qg Fock state and, consequently, is a cubic functional of the collective nuclear glue. 

The above three components of the dijet cross section differ by more than the degree of 
the nonlinearity. The coherent diffractive mechanism and the excitation of quark-gluon 
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dijets in the same color representation as the incident quark are explicitly calculable in 
terms of the collective nuclear glue of Eq. (fH4*|) which is defined for the whole nucleus. 
This is not the case for the excitation of quark-gluon dijets in higher color multiplets. It 
is proportional to the unintegrated gluon density in the free nucleon. The coherent initial 
state interaction, before the excitation of higher color multiplets at the depth /3 of the 
nucleus, must be described in terms of the unintegrated collective glue (|99|) defined for 
the slice (3 of the nucleus. Coherent distortions of the qg wave function are complemented 
by incoherent broadening of the incident quark transverse momentum distribution in the 
same slice of the nucleus. Likewise, the final state interactions after the excitation of 
higher multiplets must be described in terms of the unintegrated collective glue defined 
for the slice (1 — (3) of the nucleus. This reinforces the point |4] that hard processes in a 
nuclear environment can not be described in terms of a nuclear gluon density defined for 
the whole nucleus, as it was advocated, for instance, within the Color Glass Condensate 
approach Q. Further, hes.des the collective nuclear gl ue denned for color-slnglet 
quark-antiquark dipole, there emerges a new nuclear gluon density which depends on the 
Casimir operators of higher quark-gluon color representations, i.e., gluon field of the nu- 
cleus must be described by a density matrix in the space of color representations. Based 
on a comparison of the excitation of quark-gluon dijets in quark-nucleon collisions to the 
excitation of quark-antiquark dijets in DIS and gluon-nucleus collisions, we formulated 
four universality classes for nonlinear fc^-factorization. 

The representation for the dijet cross section similar to our master formula (|14|) has 



been discussed recently by several authors 



but these works stopped short 



of the solution of the coupled-channel intranuclear evolution for the for 4-parton state. 
Although major ingredients for the diaepnalization of the four-body S-matrix are found 
in our earlier work on dijets in DIS j4|, the case of the qg dijets has its own tricky 
points. For this reason, we felt it imperative to present full technical details of this 
diagonalization. 

The emphasis of the present communication was on the formalism, the numerical 
applications will be reported elsewhere. The nuclear coherency condition, x < xa ~ 
0.1 ■ v4~ 1//3 , restricts the applicability domain of our formalism to the forward part of the 
proton hemisphere of pA collisions at RHIC. Although the required coherency condition 
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does not hold for the mid-rapidity dijets studied so far at RHIC |38[, our predictions 
could be tested after the detectors at RHIC II will be upgraded to cover the proton 
fragmentation region LtI ]. 

This work was partly supported by the grant DFG 436 RUS 17/101/04. 
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